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ABSTRACT 

We derive the criteria for deflagration to detonation transition (DDT) in 
a Type la supernova. Tlie tlieory is based on tlie two major assumptions: (i) 
detonation is triggered via the Zeldovich gradient mechanism inside a region 
of mixed fuel and products, (ii) the mixed region is produced by a turbulent 
mixing of fuel and products either inside an active deflagration front or during 
the global expansion and subsequent contraction of an exploding white dwarf. 

We determine the critical size of the mixed region required to initiate a 
detonation in a degenerate carbon-oxygen mixture. This critical length is much 
larger than the width of the reaction front of a Chapman- Jouguet detonation. 
However, at densities greater than ~ 5 x 10^ g cm~^ , it is much smaller than 
the size of a white dwarf. We derive the critical turbulent intensity required to 
create the mixed region inside an active deflagration front in which a detonation 
can form. We conclude that the density ptr at which a detonation can form in 
a carbon-oxygen white dwarf is low, less than 2 — 5 x 10^ g cm~^ , but greater 
than 5 x 10^ g cm~^ . 

Subject headings: detonation, deflagration, supernovae, combustion 
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1. Introduction 

The standard model for Type la supernovae (SN la) is an explosion by a process of 
thermonuclear combustion. One-dimensional calculations show that the two possible modes 
of this combustion are either supersonic detonation or subsonic deflagration ( [Arnett 1969| ; 



Hansen & Wheeler 1969; [Ivanova et al. 1974| ; |Nomoto et al. 1976| ; |Nomoto et al. 1984|) . 



There is evidence, both theoretical and empirical, that SN la involve a transition from a 
deflagration to a detonation, or DDT (Khokhlov 1991a,b; [Arnett fc Livne 1994|; Hoflich, 



Khokhlov & Wheeler 1995; |H5flich 1995| ; [Wheeler et al. 1995| ; [H5flich fc Khokhlov 1996[ 



and references therein). Pure detonation models tend to produce an excess of iron peak 
elements, rather than the intermediate mass elements that are observed in the outer layers 
of SN la. Models in which DDT occurs can account for the properties of SN la, including 
the observed correlated dispersion in light curve amplitude and decay rate ([Hofiich et al 



1995[ ; [Hofiich fc Khokhlov 19961) . In such models, the burning front starts as a deflagration 



that makes a transition to a detonation directly (the standard "delayed detonation" 
model), or after the star has expanded due to partial burning and then recompressed (the 
"pulsation delayed detonation"). In these models, the density at which DDT occurs is a 
free parameter. The comparison of these models with observations suggests that DDT must 
occur at relatively low densities, of the order of 10^ g cm~'^ . In this paper we address some 
of the key physical issues of DDT. In particular, we address the physical mechanism by 
which DDT may occur in supernovae, and why it might occur at low densities. The process 
of DDT is not well understood in either supernovae or in terrestrial combustion. Many of 
the principles presented here apply to terrestrial conditions as well. 

The analysis we present is based on a theory of unconfined DDT (Khokhlov, Oran 
& Wheeler 1996b). It is assumed that the inherent mechanism of DDT is the Zeldovich 
gradient mechanism ([Zeldovich et al. 1970| ) in which ignition in a region of induction time 
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gradient leads to the development of a reaction zone and a shock that coalesce to form a 
detonation. In this theory, it is presumed that a mixed region with gradients of temperature 
and composition, a combination of which results in an induction time gradient, is created 
within a deflagration region by turbulence that microscopically mixes hot products and cold 
fuel. Two quantities must be determined: the size of the critical region, Lc, required for a 
self-propagating detonation to develop, and the amount of mixing (the intensity or velocity 
of turbulence) required to produce that critical length of mixing. 

In supernovae, there arc two situations in which DDT could occur. First, DDT could 
occur directly, when turbulent deflagration alone creates a large mixed region which ignites. 
For this process to occur, the flame surface associated with laminar burning has to be 
destroyed by turbulence. When the laminar flame front is destroyed, reactions will still 
continue at rates depending on the density and temperature, but microscopic chemical and 
thermal mixing of fuel and product occurs. The second way to initiate DDT is for the 
deflagration to die due to expansion before the DDT occurs. If insufficient energy has been 
liberated, the star will expand and then contract. Mixing of fuel and products, and the 
subsequent DDT may occur during the contraction phase. 

In this paper, we first examine the mixing process and determine the thermodynamic 
properties of the mixture as a function of the degree of mixing. We show that these 
thermodynamic properties can be characterized by a single parameter (e.g., the mass 
fraction of fuel) and that the history of mixing, that is, whether mixing occurs before, after, 
or during a pulsation phase does not change the thermodynamic properties of the mixture 
significantly. We derive the critical length scale to initiate a detonation by the Zeldovich 
mechanism. Next, we evaluate the turbulent intensity necessary to break a laminar fiame 
and to create the conditions for DDT. This required the calculation of the laminar flame 
speed and laminar flame thickness over a wide range of densities, 3 x 10^ — 10^° g cm~^ . 
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Finally, we estimate the critical density for DDT based on the existing three-dimensional 
simulations of turbulent deflagration in Type la supernovae, and discuss implications of the 
results for supernova models. 



2. Thermodynamics of mixing 

Consider a subsonic mixing process that occurs at approximately constant pressure, 
and assume that no reactions occur during mixing. The conservation of mass and enthalpy 
during mixing gives the thermodynamic parameters of the mixture, 

KTrn, Pm) = Xf h{Tf, Pf) + (1 - Xf) h{Tj„ p^) , (l) 

1 Xf il-Xf) 

— = ^ + ^ ^ , (2) 

Pm Pf Pp 

where h is the enthalpy; subscripts f, p, and m stand for fuel, product, and mixture of 
products and fuel, respectively; Xf is the fuel mass fraction in the mixture. We assume 
that the state of the products, p, is determined by conservation of enthalpy during constant 
pressure burning in a deflagration wave, 

MTp,Pp) = Mr/,P/)+Ag, (3) 

P(Tp,pp)=P(r;,p;)= const, (4) 
where Ag is the energy released by nuclear reactions. 

The value of Ag depends on the density at which the laminar flame is either quenched 
by turbulence or quenched by expansion, and is not known exactly. If the final stage of 
deflagration takes place at very low densities, say, less than lO^g cm~^, only carbon would 
be able to react in a deflagration wave. Temperatures resulting from carbon burning would 
be so low that oxygen would not react. In this case, we expect Ag ~ 3.7 x 10^^ ergs g~^, 
and the products of burning would consist predominantly of half silicon and half oxygen by 



-6- 



mass. At higher densities, both carbon and oxygen will burn to form silicon group elements. 
For the latter case, we pick a typical value of Ag ~ 6 x 10^^ ergs g~^. Below we refer to 
these two choices as Case A and B, respectively. 

Now consider mixing and changes in pressure due to global adiabatic expansion or 
contraction. For a perfect gas, the order of mixing and expansion or contraction does not 
change the final state of mixed material. For the equation of state of white dwarf material, 
this is not an exact statement, but we show that it is a good approximation. 

Figure 1 shows the results of the following numerical experiment using the degenerate 
matter equation of state described in the next section. A half carbon, half oxygen fuel at 
density 5 x 10® g cm~^ (point / in Figure 1) was burned at constant pressure to a half silicon 
and half oxygen mixture (to reach point p) and then mixed with Xf — 0.5 at constant 
pressure to yield the conditions at point nio- Alternatively, both fuel and products were 
expanded adiabatically by a factor of 100 (to points /', p'), and then mixed at constant 
pressure to produce the conditions at point m'. Finally, the mixture was returned along the 
adiabat to the initial pressure. The final state, point m, is very close to nio, even under 
this extreme adiabatic expansion and contraction. The difference in temperature between 
points m and m' is 5%, and the difference in density is 7%. We conclude that, to a very 
good approximation, the thermodynamic properties of a mixture do not depend on the 
history of the mixing and expansion processes. We make that assumption in the following 
analysis and neglect the order in which mixing and expansion or contraction occur. 

Mixing cold fuel and hot products raises the entropy of the mixture, so that the 
mixture will ignite at lower density than the fuel alone. For cases A and B, Figure 2 gives 
the density at which the energy-release rate in the fuel-product mixture is one-tenth of the 
maximum energy-release rate for isochoric burning at this density. The maximum release 
rate occurs after ignition, as the temperature rises and before fuel is consumed. Once the 
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mixture is at this density, the induction time is very short, only ten times the characteristic 
time required to burn all of the fuel. The figure shows that for a pure fuel to ignite, it 
needs to be compressed to very high densities. However, even for fuel mixed with rather 
small amounts of hot products, the ignition density decreases drastically, and is ~ 10^ — 10^ 
g cm~^ for a wide range in fuel fraction. If matter is returning to high densities after a 
pulsation, a mixture containing a noticeable amount of products cannot be compressed to 
densities significantly in excess of 10'' g cm~^ before it ignites. For the pulsation delayed 
detonation models, therefore, the thermodynamic considerations place a limit on the 
conditions at which DDT can occur: DDT can only occur at densities near or below 10^ g 
cm~^. For the standard delayed detonation models, however, no such limit can be placed by 
thermodynamic considerations. For both cases, the kinetics and hydrodynamics of ignition 
need to be considered to determine whether DDT can occur at thermodynamically allowed 
conditions. 



3. Critical length scale for DDT 

We now determine the critical size of the mixed region capable of triggering a 
detonation. It has been shown ( [Khokhlov et al. 1996b| ) that there are two processes that 



are necessary for DDT to occur by the Zeldovich mechanism. First, a spontaneous wave 
must be generated and a shock-reaction complex must be formed. This complex forms 
somewhere in the middle of the mixed region where there is both enough fuel and high 
enough temperature. Second, this complex must be able to survive the propagation into 
the unburned material. 



To estimate the critical size of this region, Lc, we consider, as in ([Khokhlov et al. 



1996b|) , a nonuniform region created by mixing the products of isobaric burning and fresh 
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fuel, such that there is a hnear spatial distribution of fuel fraction Xj, 

f x/L, < X < L 
Xf{x) = , (5) 

[1, x> L 

where L is the size of the mixed region. This creates a region of oppositely directed 
gradients in temperature and concentration of reactants. Initially, the velocity of the 
material is zero, and the pressure Pq is constant everywhere. The boundary conditions at 
X = are reflecting walls (symmetry conditions). 

The system is described by the Euler equations coupled with the nuclear kinetic 
equations 

1^ + V ■ (pU) = , (6) 
^ + V-(pUU) + VP = 0, (7) 
+ V-((P + P)U) =0, (8) 

^ + (U-V)3^ = 7^, (9) 

where the total energy density E = Et + pU^ /2 + pNaQ ■ y is the sum of the 
thermal, kinetic, and nuclear energies, U is the fluid velocity, Et is the thermal energy 
density, 3^ is the vector in the composition space of mole fractions of nuclei, Q is the 
vector of corresponding binding energies per nucleon, and 7^(p, Et, y) is the vector 
of net rates of change of nuclear species with time; Na is the Avogadro number. 
To describe thermonuclear burning, an alpha-nuclei reaction network is used with 
3^ = {^He, i^C, 1^0, 20^6, 24m^, 28^2,32 s, ^Mr, 40 Ca, ^^Ti, ^^Cr, ^^Fe, ^'^Ni). The equation 
of state includes contributions from Fermi gases of electrons and positrons, radiation, and 
ions. Further details of the input physics can be found in ( [Khokhlov 1993| ). 

The system of equations is integrated numerically using a one- dimensional version 
of the time-dependent, compressible fluid code ALLA based on the Piecewise Parabolic 
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Method ( ICollela fc Woodward 19"8l ; [CoUela fc Glaz 1985 ), implemented as lagrangian 



step plus remap. Details of the implementation are given in ( [Khokhlov 1995| ; [Khokhlov 



|et al. 1996a|) . The nuclear reactions are coupled to fluid dynamics by time-step splitting. 



The kinetic equations are integrated together with the equation of energy conservation 
using a stiff solver with adjustable substeps to keep the accuracy better than 1%. Most 
computations were done using a grid with 1024 equidistant computational cells. The 
convergence of the numerical solutions was tested for selected cases by varying the number 
of cells from 1024 to 4096. Computations were carried out in both planar and spherical 
geometries. 

Computations were performed assuming that the products of burning that formed 
the mixture (state p in the thermodynamic relations given in equations 1 to 4) are 
characterized by the value of the energy released by the deflagration Aq = 4 x 10^"^ ergs 
g~^, or Aq = 3 X 10^"^ ergs g~^. These values are close to Case A of § 2 and have been 
selected to test the sensitivity of the results to Aq. For Aq = 4 x 10^"^, the composition of 
products was taken to be pure "^^Si . For the second value, the composition of products was 
taken to be half "^^Si and half ^^O by mass. The system was prepared in an initial state 
and then allowed to evolve in time, flrst until ignition takes place, then until the formation 
(or failure) of detonation, and flnally to the time when the generated detonation or shock 
leaves the computational domain. The results of these computations are not sensitive to 
the assumed energy release and composition. They mainly depend on the density at which 
explosion occurs, and on the size of the explosion region. 

Figure 3 is an example of a computation in which a detonation successfully forms. 
The initial density of the fuel po = 1-0 x 10^ g cm~'^ and the size of the mixed region is 
L = 5.0 X 10^ cm. Figures 3a through 3e show proflles of density, temperature, pressure, 
velocity, and carbon mole fraction, respectively, for selected times. In the initial state, the 
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density and carbon mole fraction of the mixture increase, and the temperature decreases as 
a function of position x. Because of such opposing temperature and concentration gradients, 
the ignition does not occur first at the origin, x = 0, but is shifted to some point x > Q 
where the energy generation rate is maximum. Prom this point, the ignition spreads out 
supersonically as a spontaneous wave. A detailed discussion of spontaneous burning is given 
in Khokhlov et al. (1996b). The overpressure and material velocity in the spontaneous wave 
grow with time as the phase velocity of the burning front decreases. When the speed of the 
spontaneous wave approaches the speed of sound, a shock wave emerges, and a shock wave 
- reaction complex forms. The profiles 1 to 3 illustrate the spontaneous wave propagation. 
There is no shock present at this stage. Matter is continuously compressed and accelerated 
in the process of burning, and expands and slows down afterwards. The shock - reaction 
complex forms between the profiles 3 and 4. The complex continues to propagate through 
the matter in which density increases and temperature decreases (profiles 4 and 5), but the 
strength of the complex grows. Finally, the complex passes into the fully unburnt cold fuel 
in the form of a slightly overdriven detonation (profile 6), and begins its relaxation to a 
Chapman- Jouguet state. 

Figures 3f and 3g show the nuclear composition for the times corresponding to profiles 
2 and 7, respectively. These figures illustrate that at densities ^ lO^g cm~^, the important 
reactions are ^^C + ^^C — > ^°iVe , "^^Mg , "^He , and the reactions on ^°iVe and "^^Mg that 
transform the products of ^^C + ^^C burning into the "^^Si -group elements. Reactions 
involving oxygen are less important. The initial amount of oxygen in the mixture remains 
practically intact. The energy release is approximately Ag ~ (3.5 — 4) x 10^^ ergs g^^, 
consistent with the value of the energy release for the products of defiagration adopted 
above, and with the assumptions about the energy release made in § 2. Figures 4a and 
4b show a successful initiation of a detonation at a density five times higher, po = 5 x 10^ 
g cm~^ . At this density, burning leads to higher temperatures due to a decrease in the 
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specific heat of degenerate matter. As a result, oxygen burns to silicon rather quickly 
(Fig. 4b), and the energy of the explosion, Aq ~ 7 x 10^7 ergs g~^, is almost two times larger 
than in the previous case po = lO^g cm~^. 

The shock-reaction complex can survive propagation down the temperature gradient 
and can grow into a fully developed detonation only if conditions in front of the complex 
change slowly enough, that is, if L is large enough. If the variations of preshock density and 
temperature with distance are too fast, the reaction zone structure is unable to adjust to 
the decreasing reaction rates behind the shock, and is unable to compensate for the energy 
required to shock denser fuel. 

Figures 5a to 5f illustrate the failure of the initiation of a detonation at po = lO^g cm~^ 
when the size of the mixed region was 20 percent smaller, L — 4 x 10^ cm, than in the 
case presented in Fig. 3. Here, too, the spontaneous wave, and the shock - reaction 
complex form (profiles 1 to 3). However, the complex does not survive the propagation 
down the temperature gradient. The shock and reaction separate. The shock passes into 
the unburned matter leaving the reaction wave behind. It is the condition that the complex 
must survive the propagation down the temperature gradient, not the formation of the 
complex alone, that determines the critical size of the mixed region. 

Figure 6 summarizes the results of a series of many computations performed to 
determine Lc, the critical size of the region capable of triggering a detonation, as a function 
of fuel density. The value of Lc is a very sensitive function of density. It is virtually 
independent of assumptions made about the energy release Aq and composition of the 
products formed during the preceeding deflagration stage. The critical length depends, 
however, on the energy release during the explosion itself. As was mentioned earlier, the 
energy release increases with increasing background density, which makes the dependence of 
Lc on density very strong. The values of Lc obtained in spherical and planar geometries are 
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practically the same, because the radius of curvature of a spherical front is approximately 
the same as the critical length which is much larger than the detonation wave thickness. 

The simulations show that the spontaneous wave grows to an appreciable strength, and 
a shock - reaction complex forms in a part of the mixed region with a high fuel fraction, 
Xf ~ 0.8 — 0.9. Thus, it is this mixture that is important for DDT. As shown in Figure 2, 
explosion of such a mixture can take place at densities < 10 ''g cm~^ (Fig- 2). On the other 
hand, as shown in Figure 6, the critical length for DDT is a sensitive function of density; 
it increases steeply at low densities and it is of order 3 x 10^ — 4 x 10^ cm at densities of 
5 X 10^ — lO'' g cm~^. This length is a factor of 10^ — 10^ larger than the size of a detonation 
reaction zone. However, even at these low densities, this length is small compared to the 
radius of a white dwarf. The length scale L^. is small enough to be formed by turbulent 
mixing in a relatively small portion of the supernova ejecta for p ^ 5 x 10^ g cm~^. 

4. Conditions on turbulent intensity 

Here we discuss the conditions necessary to produce a mixed region of size Lc for both 
the direct and the pulsation modes of delayed detonation. First consider DDT in the direct 
mode. A fundamental assumption is that turbulence must destroy the flame front so that 
burned and unburned materials can mix microscopically. To break the flame front, the 
turbulence velocity Ub, at a scale comparable to the thickness of the laminar front, must 
exceed the normal speed of the laminar flame, 

Ub^K 5,, (10) 

where Se is the laminar flame speed, Se is the laminar flame thickness, and i^' ~ 1 is 
a coefficient that describes the ability of the laminar flame to survive in a turbulent 
environment ( [Khokhlov et al. 1996b|) . 
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The mechanism for breaking the flame is flame stretching, which includes the effects of 
strain and curvature. These effects depend on the properties of the particular flame and 
on the spectrum of the turbulence. For terrestrial flames, the coefficient K significantly 
exceeds unity according to both theory and experiments (Poinsot, Candel fc Trouve 1996| ; 



Roberts et al. 1993|) . For the thermonuclear fiame, the value of K is not known exactly. 
The analysis of the nonlinear stabilization of the fiames in a turbulent field suggests = 8 
( [Khokhlov 19951 ). The numerical simulations of the fiame presented in ( [Khokhlov 1995| ) are 



consistent with a large value of K. For the subsequent analysis in this paper, we use the 
two values K = 1 and K = 

Figure 7 shows the laminar fiame speed Si and thickness 5i for a 0.5C + 0.5O mixture 
computed for the range of densities 3 x 10^ — 10^° g cm~^ , as described in the appendix. 
In the same figure, we show Si according to Timmes and Woosley (1992) who computed 
the laminar fiame speed at high densities, p > 2 x 10^ g cm~'^ . Our results are in good 
agreement with theirs at these densities. At low densities the fiame speed decreases with 
density more rapidly than at high densities, and reaches a very low value of 5^ ^ 3 x 10^ 
cm/s at p ~ 3 X 10'' g cm~^ . The computed values of 5*^ give Ub assuming K = 1. The 
curve = 8Si shows the required turbulent intensity ii K = 8. 

To relate Ub to the intensity of turbulent motions on larger scales, additional 
assumptions are required about the spectrum of turbulent motions inside the white dwarf 
during the explosion. The most favorable assumption for DDT is that the turbulent 
spectrum is Kolmogorov. By assuming that a Kolmogorov spectrum is established, we 
ignore the time delay required to establish the steady-state cascade of turbulent energy, and 
ignore freeze-out of turbulence due to the expansion of the star. Under the assumption of a 
Kolmogorov spectrum, the speed of turbulent motions at the scale of the laminar front 6i is 
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related to the speed of turbulent motions on larger scales L > Si as 

C.,=%(f)'". (11) 

During the deflagration phase of the explosion, large portions of the star will be 
occupied by the very convoluted flame. To be speciflc, we pick a typical value L = L/ ~ 10^ 
cm for the maximum size of the turbulent flame region ( [Khokhlov 1995| ). The turbulent 
region may be somewhat smaller in the beginning of the explosion, and somewhat larger 
at the end, but this value is of the correct order of magnitude. If we assume that the 
deflagration is strong enough to unbind the star, we have to assume that Lj is large 
enough. The large-scale turbulence in a supernova is driven by buoyancy. If we neglect 
the effects of freeze-out due to expansion of the star, we expect Ulj to be of order of the 



Rayleigh- Taylor characteristic velocity at this scale Ulj, ~ tZ/j-T — 0.5y where g is 
the effective gravitational acceleration. In a star in hydrostatic equilibrium, or close to it. 



the value of JgLj does not exceed the speed of sound. In equilibrium, P/R ^ pg, where 



R is of order the stellar radius. From this it follows that ~ J P/ p ^ VgR and hence 



gLf < as for Lf < R. For our estimates we take the value of = Ur-t = 0.5as. The 
values of the turbulent velocity at the scale of the flame front estimated from equation (11) 
using these values is plotted in Figure 7. The intersection of the curve Us^ with the curves 
Si and gives estimates of the transition density ptr < (4 — 8) x 10^ g cm^^ below which 
DDT may occur. We emphasize that this estimate of ptr is made neglecting any time delay 
in establishing the Kolmogorov spectrum, and neglecting freeze-out of turbulence due to 
expansion of the star. 

In order to take the effect of freeze-out into account, three-dimensional simulations of 
the explosion of the entire star are required. To date, only one such simulation has been 
performed ( Khokhlov 1995| ), and the analysis presented in that paper suggests that the 



expansion freezes scales exceeding approximately L ^ 10^ — 10'' cm, and that the turbulent 
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burning velocity at these scales does not exceed Ul ^ 10^ — 10'' cm/s. If we take from this 
range the most favorable values for DDT, Ul = 10^ cm/s and L = 10^ cm, equation (11) 
gives another estimate of the turbulent velocity at the scale of the flame front, [/|^, which is 
also shown in Figure 7. The intersection of the curve f/|^ with the curves and Si gives, 
we believe, a more realistic estimate of the transition density ptr ^ (2 — 5) x 10'^ g cm~^ . 

There are other effects that may influence the transition density. Self-turbulization 
of the flame on scales less than Lf might be possible, as is preconditioning of a mixture 
by shocks generated by the turbulence itself. The latter effects might somewhat facilitate 
DDT, but not to a great extent (Khokhlov et al. 1996b). Self-turbulization due to the 
Landau-Darrieus instability is limited by nonlinear effects. Shock preconditioning requires 
a very strong shock to be effective. We believe that there are two major unknown factors 
in estimating the value of ptr^ the spectrum of turbulence in the exploding star, and the 
extent to which the flame can survive turbulent stretch on small scales (coefficient K). 

Now consider the case in which the conditions for DDT do not occur in the deflagration 
phase, and the deflagration is extinguished by the expansion. In this case, burning is 
quenched by expansion and the products and fuel mix freely. There is no laminar flame 
front and there are no nuclear reactions. This situation is still Rayleigh- Taylor unstable. 
Previous estimates ( Khokhlov 1991a ; Arnett fc Livne 1994 ) have shown that the mixed 



R-T region may be on the order of 10^ — 10*^ cm, which is much larger than at 
p ^ 5 X 10^ — lO'^g cm~^. With a Kolmogorov cascade, the eddy turnover time scales as 
A^/^, where A is the size of the eddy. When the star is compressed, the turbulent motions 
on all scales will be enhanced and mixing should proceed vigorously (this is opposite to 
the effect of freeze-out expected during the expansion of the star). It is thus reasonable to 
expect microscopic mixing to be complete before the contraction phase.Q This suggests that 



^We have invoked microscopic mixing in this model. There is considerable evidence that 
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the pulsation mechanism can lead rather naturally to conditions where DDT is unavoidable 
in the recompression phase. 

5. Conclusions 

We applied a theory of DDT in unconfined conditions ( jKhokhlov et al. 1996a| ) to DDT 
in Type la supernovae. The two basic assumptions of this theory are: 1) the gradient 
mechanism is the inherent mechanism that leads to DDT in unconfined conditions, and 
2) the mechanism for preparing the gradient of induction time is turbulent mixing that 
requires breaking or quenching the flame front. 

Using a series of numerical simulations, we determined the minimum critical size of 
the mixed region, Lc, required for DDT (§ 3). We find that Lc, though much larger than 
the detonation wave thickness, is at least two to three orders of magnitude smaller than 
a white dwarf radius for p > 5 x 10^ — 10^ g cm~^ . Thus, a mixed region of critical size 
required for triggering a detonation can, in principle, be formed by turbulent mixing in a 
small portion of a supernova. 

There are two possible modes of detonation formation in a Type la supernova: a direct 

mixing occurs in supernova ejecta but is not always microscopic. Several lines of evidence 
strongly suggest that the mixing in SN 1987A was macroscopic, but not microscopic. We 
note, however, that the physical conditions here and in a collapse environment like SN 1987A 
are different. The instability leading to mixing in SN 1987A is expected to be a Richtmeyer- 
Meshkov instability in which the density /pressure inversion is only temporarily induced at 
composition/density boundaries due to shock passage. In the thermonuclear explosions that 
we consider here, the Rayleigh- Taylor instability is intrinsic and long lasting. There should 
be ample time to mix to small length scales. 
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mode when DDT occurs inside an active deflagration front, and a pulsation mode which is 
possible if DDT does not occur directly, and if the energy released during the deflagration 
phase is not enough to unbind the star. Then the star will experience a global pulsation, 
and mixing will take place between the cold fuel and ashes of the dead deflagration. The 
requirements on the intensity of turbulent mixing needed to form a mixed region of critical 
size Lc are different in these two cases. 

For a direct mode of DDT (delayed detonation model), the intensity of turbulence on 
small scales of the order of the laminar flame thickness must be high enough to stretch 
and break the surface of the laminar flame front (§4, equation 10). We computed the 
laminar flame velocity and thickness in a wide range of densities 3 x 10^ - 10^° g cm , 
and estimated the critical intensity of turbulence required to break the flame at the scale of 
a laminar flame thickness. 

We also tried to estimate the value of the transition density ptr at which DDT should 
occur during a Type la supernova explosion. Using available information about the 
intensity of turbulence inside the exploding star, we conclude that a realistic estimate of the 
transition density for the direct delayed detonation explosion is ptr ^ (2 — 5) x lO^g cm~^. 

However, even under the most favorable and unrealistic assumptions (§ 4), ptr does not 
exceed 10^ g cm~^ . 

In the pulsation mode (pulsating delayed detonation model), microscopic mixing 
throughout large areas can occur because there is no flame front separating burned 
and unburned matter, and, thus, fuel and products can mix freely. We considered the 
thermodynamics of the mixing process in a supernova (§ 2) to demonstrate that the order 
of mixing and global expansion or contraction of the supernova is not important. We also 
found that even a small fraction (^ 1%) of high-entropy products of burning mixed with 
cold, low-entropy fuel increases the entropy of the mixture such that it will ignite at very 
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low densities. We conclude that in this mode, the portions of the mixture of size comparable 
to or exceeding Lc would unavoidably detonate around ptr — 10^ g cm~^. 

Despite various details of DDT that require future consideration, our main conclusion 
is that DDT in Type la supernovae should occur in the density range 5 x 10^ to a few 
times 10^. At higher densities, 10^ — 10^ g cm~^ , either the required turbulent velocity 
cannot be reached (direct mode of the delayed detonation) or the material inevitably ignites 
before such high densities are reached (pulsating mode of the delayed detonation). At 
lower densities, the size of the critical region becomes too large. This conclusion agrees 
with the results of global modeling of Type la supernovae, where the transition density, ptr, 
at which DDT occurs is a free parameter. Global modeling suggests that, in order to fit 
observed hght curves and spectra of SNIa, the transition density should be ptr — 10'' g cm~'^ 
( IKhokhlov 199Ta| ; [Woosley fc Weaver 1994] ; |Hotlich et al. 1995| ; |Nomoto et al. 1995| ; [Hotlich 



[T995t |H5tlich fc Khokhlov 19961) . 
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A. Properties of laminar flames in degenerate C-O mixtures 

The velocity of a laminar flame Si is computed as the eigenvalue of the diflerential 
equations describing the structure of a steady-state burning wave ( [Zeldovich et al. 1985|) . 
We assume that the wave propagates steadily and that burning occurs at constant pressure 
since the flame speed is much less than the sound speed. The parameters of the fuel and 
the products of burning are designated below by subscripts and 1, respectively. The mass 
flux is constant through a steady wave, = pu = poSi = piUi. The equation of energy 
balance inside the wave is 



where cp is the speciflc heat at constant pressure, x is the thermal conductivity, Na is the 
Avogadro number, Q is the vector of binding energies of nuclei, TZ is the vector of net 
reaction rates (Section 3), and x is the distance across the wave. The terms on the right 
side of (A.l) describe the temperature changes due to the heat conduction and nuclear 
reactions. Nuclear mole fractions y are changing inside the wave due to nuclear reactions, 



Ion diflusion is ignored in (A. 2) since, in a degenerate gas, the mean free paths of ions are 
much smaller than the mean free paths of electrons and photons. The density p(T, 3^) inside 
the wave is related to T and y via the constant pressure condition 




(Al) 



ax 



{A.2) 



P{p,T,y) = Po. 



{A.3) 



The boundary conditions for Eqs.(A.l) to (A.3) are 



at a; = — oo : T = Tq, p = po, y 




(AA) 



and 



at X = oo : T = Ti, p = pi, y 




. 



(A.5) 
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The points x = ±00 are singular points of equations (A.1)-(A.3). The integral curve which 
passes through these singular points represents the solution for a steady flame front. We 
designate the value of Jm corresponding to this integral curve as J*. 

For low temperatures, T < T^, where is some value less than 2 x lO^K, we can 
neglect the energy generation term in (A.l) in comparison with the term describing the 
heat transport. At these temperatures, (A.l) can be rewritten as 

dW _ JmCp 

where W = X%:- The corresponding initial condition for (A. 7) is W{Tq) = 0. The 
integration of (A. 7) from Tq to gives a certain value of = W{T^), from which we find 
the corresponding value of (^)^ • From this point, we integrate the full equation (A.l) 
with the energy generation term included. 



The value J* is found using the method of trial and error (|Zeldovich et al. 1985|) 



The integral curves for 7^ J* diverge from the J* curve when x ^ 00. For < J* 
the asymptotic (x 00) behavior is T ^ —00, ^ —00. For > J* the asymptotic 
behavior is T — > 00, ^ 00. We use the difference in the behavior of the integral 
curves for discriminating between the cases Jm < J* and Jm > J*- The next trial value 
of Jm is found as the average of the two most recent trial values lying on different sides 
of J*. It requires ~ 20 iterations to obtain the value of J* with an accuracy of 10"'^. 
Then the laminar flame speed is calculated as Si = J*/po- The values of 61 are estimated 
approximately as the thickness of a zone where most of the nuclear energy is released and 
the CO mixture is converted into Si-peak elements. Computations are done for the initial 
mixture of equal masses of ^^C and ^^O with a small amount of "^"^Ne added in order to get 
a neutron excess 77 = 0.002. The initial temperature is Tq = lO^K, and we used = 10^ 
K. Flame parameters are not sensitive to variations of Tq and if these temperatures are 
below 2 X lO^K. 
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In computations of Si and Si, the following 133 species reaction network is used: n, 

p, ^He, 13 AT, 160, 20,22 JVe^ 23 jy^^ 23-26 27^;^ 27-32^^^ 30-33p^ 31-36^^ 35-37^/^ 

36-41 Ar, 39-43;^^ 40-460^^ 43-47^^^ 44-50y,^ 47-52^^ 48-56^^^ 51-60 52-62^^^ SS-ei^o, 

56-64 jY^^ 57-65(jy^^ 59-66 ^'j^ rpj-^g forward reaction rates are taken from Thielemann's reaction 
rate library ( P?hielemann, Arnould fc Truran 1987|) . The binding energies and partition 
functions were taken according to Woosley et al. (1987). The backward reaction rates are 
calculated according to the principle of detailed balance. Screening corrections are applied 
to the most important reactions, ^^C +^2 c, ^^C O, ^^O O and 3 '^He ^^2 o, which 
control the thermonuclear energy release at low temperatures. Even for these reactions 
screening factors were found to be relatively small {E <^2 — 5) for conditions where nuclear 
reactions significantly contribute to the right side of Eq.(2). Most of the nuclear reactions 
occur at high temperatures and the corresponding screening factors ^ 1.5 are small 
in comparison with the expected uncertainties in the nuclear reaction rates (~ 2 — 3). 
Screening corrections are taken from Yakovlev and Shalybkov (1989). 

The total thermal conductivity is calculated as the sum of electron and photon 
conductivities 

X = Xe + Xf, {A.8) 



with x-y = 4:acT^ /Sup. The expression for the photon opacity ( [Paczynski 1983|) 



4 r / / T \0-86\ 2 



(1 + 2.7 10iip/T2) \ V4.5 108 
approximates the low density (^ 6 lO^g cm^^) results of Buchler and Yueh (1976) that take 
the degeneracy of electrons and relativistic scattering into account. Ye is the electron mole 
fraction. At higher densities, the contribution of the photon opacity to the total opacity is 
very small. Our results at high densities are in a good agreement with those of Timmes & 
Woosley (1992) who used different approximations for the photon thermal conductivity. 

The electron thermal conductivity can be expressed in terms of the effective electron 
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coUisional frequency, u^i ^ (Yakovlev and Urpin 1980) 

Xe = 4.09 10^ T f { i^") ^ {A.m 

where x — p-p/nieC ~ 1.009 10~^(pFe)^^^, Pf is the electron Fermi momentum, and is a 
sum of ion-electron and electron-electron coUisional frequencies. 

The electron - ion coUisional frequency is calculated using the mean ion approximation 
following Urpin & Yakovlev (1980) 

V,, = 1.78 lO^^yrT^ — ^, {A.V1) 

^ e 

where Y — and Z — Yf./Y are the mean ion mole fraction and the mean atomic 
number, respectively, and A is the Coulomb logarithm which can be expressed for gaseous 
and liquid ionic states as the sum of the logarithm obtained in the Born approximation and 
the non-Born correction (Yakovlev and Urpin 1980, Yakovlev 1987) 



A = Ab + (A. 13) 



with 



Ae^.((M)-/^)-^, 

r = 2.275 10^ -jr {pYef^ {A.lb) 



and 



The thermal conductivity due to electron-electron collisions of relativistic electrons in 
a degenerate electron gas was calculated in the static screening approximation by Urpin 
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and Yakovlev (1980), 



^3/2 



i/ee = 0.511 J{y) s-^ , (A. 18) 

(1 + x^) ' 



where 

(A. 19) 



T 
and 

T^e = ^ = 3.307 10^ K (A.20) 

is the electron plasma temperature. Recently Yakovlev (unpublished) recalculated the 
integral entering the general expression of v^e numerically and obtained the following 

analytical expression for J(|/) 

with a — 0.113 and h — 1.247. The largest error in the fit is ~ 6% at y ~ 2. 
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Fig. 1 - Comparison of thermodynamic properties of mixtures m and mo formed by 
mixing the same products and fuel at high and at low pressure, respectively. 

Fig. 2 - Ignition density as a function of the fuel fraction. (Sohd hne) - Aq — 3.7 x 10^^ 
ergs g~^; (Dashed line) - = 6.8 x 10^'^ ergs g~^. 

Fig. 3 - Initiation of a detonation for po — 10^ g cm~^ and L — 10^ cm. Profiles 
are shown for the initial time t = 0, and for times t = 2.86 x 10~^ (1), 7.05 x 10~^ (2), 
1.40 x 10-4 (3), 1.63 X 10-4 (4), 1.90 x 10-^ (5), 2.33 x 10-^ (6), 2.76 x 10-^ (7) and 
3.27 X 10"^ (8) sec. (a) Pressure, (b) Temperature, (c) Density, (d) Velocity, (e) Carbon 
mole fraction, (f) Distribution of nuclear species at the time corresponding to profile 2. (g) 
Distribution of nuclear species at the time corresponding to profile 7. 

Fig. 4 - Initiation of a detonation for — 5 x 10^ g cm"^ and L = 1.3 x 10^ cm. 
Profiles are shown for the initial time t = 0, and for the moments of time t = 1.74 x 10-^ 
(1), 7.89 X 10-s (2), 1.38 x 10-^ (3), 1.90 x 10-^ (4), 2.41 x 10-^ (5), 3.38 x 10-^ (6), 
4.31 X 10"^ (7) and 4.96 x 10"^ (8) sec. (a) Pressure, (b) Distribution of nuclear species at 
the time corresponding to profile 7. 

Fig. 5 - Failed initiation of a detonation for po = 1 x 10^ g cm"^ and L = 8 x 10^ cm. 
Profiles are shown for the initial time t = 0, and for the moments of time t = 5.64 x 10~^ 
(1), 1.12 X 10-4 (2), 1.68 X 10-4 (3), 2.22 x 10-4 (4), 2.77 x 10-4 (5) and 3.31 x 10-4 (6) 
sec. (a) Pressure, (b) Temperature, (c) Density, (d) Velocity, (e) Carbon mole fraction, (f) 
Distribution of nuclear species at the time corresponding to profile 5. 

Fig. 6 - Critical length for the initiation of a detonation as a function of the fuel 
density. (Circles) - Ag = 4 x 10^''' ergs g"^ for defiagration products. (Triangles) - 
Aq = 3 X lO^'' ergs g-^ for deflagration products (all but one triangle coincide with the 
circles) . 
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Fig. 7 - Velocities and flame thickness as a function of density. -S"; and 5i - laminar 
flame speed and thickness, computed in the appendix. Sf^ - laminar flame speed according 
to Timmes & Woosley (1992). - critical turbulent velocity assuming K = - 
estimated turbulent velocity at the scale of the flame front during the supernova explosion, 
assuming no turbulence freeze-out. U9 - same, but assuming turbulence freeze-out. 
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